Phonon-assisted decoherence in the production of polarization-entangled photons in a 

single semiconductor quantum dot 
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We theoretically investigate the production of polarization-entangled photons through the biex- 
citon cascade decay in a single semiconductor quantum dot. In the intermediate state the entangle- 
ment is encoded in the polarizations of the first emitted photon and the exciton, where the exciton 
state can be effectively "measured" by the solid state environment through the formation of a lat- 
tice distortion. We show that the resulting loss of entanglement becomes drastically enhanced if 
the phonons contributing to the lattice distortion are subject to elastic scatterings at the device 
boundaries, which might constitute a serious limitation for quantum-dot based entangled-photon 
devices. 
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Single quantum dots provide a viable source for en- 
tangled photons [1, 2, 3], as needed for quantum com- 
putation [4] and quantum cryptography [5]: a biexciton, 
consisting of two electron-hole pairs with opposite spin 
orientations, decays radiatively through two intermedi- 
ate optically active exciton states. If the exciton states 
are degenerate, the two decay paths of the cascade dif- 
fer in polarization but are indistinguishable otherwise. 
Therefore, the emitted photons are entangled in polar- 
ization [1]. Such ideal performance is usually spoiled by 
the electron- hole exchange interaction [6, 7], that splits 
the intermediate exciton states by a small amount and 
attaches a which-path information to the photon fre- 
quencies. This process deteriorates the photon entan- 
glement. Spectral filtering [2], energetic alignment of the 
exciton states by means of external magnetic [3] or elec- 
tric [8] fields, or growth optimization [9] allow in princi- 
ple to erase this information and to recover a high de- 
gree of entanglement. While all these schemes mask the 
which-path information to an outside observer, who could 
gain information about the photon polarization by simply 
measuring its frequency, it might also be the solid-state 
environment within which the quantum dot is embedded 
(e.g., lattice degrees of freedom or charging centers) that 
acts as an "observer" of the intermediate exciton state 
and thereby diminishes the entanglement degree. 

It is the purpose of this Letter to theoretically inves- 
tigate the loss of photon entanglement due to couplings 
to the solid state environment. We consider the situa- 
tion that the interaction of the exciton with the lattice 
degrees of freedom, i.e., phonons, depends slightly on the 
exciton polarization. This is expected for self-organized 
quantum dots owing to the strong strain-induced piezo- 
electric fields [9, 10] that are also responsible for the 
fine-structure splitting. We show that the biexciton cas- 
cade decay is accompanied by the formation of a polaron 
[11, 12, 13, 14, 15], i.e., a lattice distortion in the vicinity 
of the quantum dot, leading to a partial loss of coherence 
and photon entanglement. This effect becomes drasti- 



cally magnified for single-dot devices, such as mesas or 
microcavities. The phonons contributing to the polaron 
can scatter elastically at the device boundaries [16]. In 
consequence, the phonons continuously "measure" the 
exciton state and einselect [17] the system from a max- 
imally entangled state to a mixed one. We show that 
this effect leads to a significant reduction of the degree 
of photon entanglement, and might constitute a serious 
limitation for quantum-dot based entangled-photon de- 
vices. 




FIG. 1: (color online) Schematic sketch of the creation of po- 
larization entangled photons in a single semiconductor quan- 
tum dot through the cascade decay of the biexciton XX, as 
described in text. The two intermediate exciton states Xh 
and Xv have dipole moments oriented along x and y, respec- 
tively; is the quantum dot groundstate. The different panels 
report the excitonic wavefunctions and the distortion of the 
lattice in the vicinity of the quantum dot, which is assumed 
to be slightly asymmetric for states Xh and Xv- Hxx {Hx) 
and Vxx (Vx) denote the polarizations of the photons emitted 
in the biexciton (exciton) decay. We assume a well-defined 
polarization of the emitted photons, which can be achieved 
by placing the quantum dot into a microcavity [1, 2, 3], and 
equal magnitude of the dipole moments. The frequencies of 
the photons created in the biexciton and exciton decay differ 
by a few meV owing to the biexciton binding energy. 
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In our theoretical approach we consider a single quan- 
tum dot which is initially in the biexciton state |3) 
(two electron-hole pairs with opposite spin orientations 
[1, 2, 3], XX in Fig. 1). In a typical experiment, the 
system will eventually pass after photoexcitation or elec- 
trical injection of electrons and holes through this state. 
It thus suffices to consider this initial condition. The 
biexciton is optically coupled to the two exciton states |1) 
(Xh) and |2) (Xy ) which are polarized along x and y, re- 
spectively [7, 9, 10]. We presume that the fine-structure 
splitting due to the anisotropic electron-hole exchange 
interaction is corrected for by appropriate external mag- 
netic or electric fields [3, 8], such that the energies E\ 
and E2 are identical. In the ideal case, the cascade decay 



|3> -» \H1) + \V2) -> \HH) + \VV) 



(1) 



then generates a maximally entangled two-photon state 
[1, 18, 19], where for simplicity we have suppressed all 
phase and normalization factors. Under realistic con- 
ditions, however, the state of the emitted radiation is 
affected by a number of uncertainties, such as the pho- 
ton emission time (time jitter [20]) or environment cou- 
plings. In particular, if the two states 1 and 2 couple 
differently to the solid state environment, for reasons de- 
tailed below, the intermediate entangled photon-exciton 
state will gradually decay, with a given rate 7, into a 
mixed state. To quantify the resulting entanglement loss, 
we introduce the two-photon correlation function [21] 
G (2 \t,T) — (:Ixx(t)Ix(t + t):), describing the polar- 
ization correlations between the photons emitted in the 
biexciton decay at time t and the exciton decay at time 
t + r, respectively. It can be computed by means of the 
quantum regression theorem [18, 21]. In the two-photon 
subspace spanned by the basis HH, HV, VH, and VV, 
we then obtain 
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with r the radiative decay rate. The first exponential 
on the right-hand side accounts for the gradual buildup 
of the two-photon state due to the radiative biexciton 
and exciton decay, while the matrix describes the (condi- 
tional) two-photon density matrix. For the specific form 
given in Eq. (2) one can easily check that the entangle- 
ment of formation or concurrence [22] , which provides 
a quantitative measure of the photon entanglement, is 
simply given by e ,T . Thus, the longer the interme- 
diate exciton interacts with the environment, the more 
the entanglement becomes diminished. The overall de- 
gree of polarization entanglement 1/(1 + 7/T) is obtained 
through an average over all times t and r [18]. 

Having established the consequences of polarization- 
dependent exciton couplings, we shall now propose a mi- 
croscopic model for such asymmetry. Our starting point 
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FIG. 2: (color online) Decay of polarization amplitude 
|{2|/o(t)|l) | in presence of phonon couplings and for three dif- 
ferent temperatures. At time zero the system is prepared in 
the superposition state |1) + |2). The inset shows for compar- 
ison the time evolution of the optically induced polarization 
j(l|p(t)j0)| (initial state |0) + |1)), that is directly accessible in 
four- wave mixing experiments [12]. The gray lines correspond 
to simulations where elastic phonon scatterings are ignored, 
the dark lines to simulations with jd = 50 fieV [16] . In the 
calculations we assume an asymmetry n = 0.1 and Gaussian 
exciton wavefunctions with a FWHM of 5 nm and 1.5 nm 
along the lateral and growth directions, respectively [23]. 



is provided by the usual independent Boson model [11, 13] 
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which describes the coupling of the biexciton and the 
excitons to the phonons. Here, H = Y^%= 1 Ej\ + 
^2 X w\a\a^, with A labeling the different acoustic phonon 
modes of energy u>\, a\ are the bosonic creation opera- 
tors, and g\ the (bi)exciton-phonon couplings. Since the 
biexciton approximately consists of two excitons with op- 
posite spin orientations, we set g^ = g\ + g\. In an 
ideal quantum dot the coupling constants g\ = g\ are 
given by the form factors of the exciton wavefunctions 
[13, 23] and the phonon coupling constants of the bulk 
semiconductor. These are entirely real for deformation 
potential coupling and entirely imaginary for piezoelec- 
tric couplings [11]. For self-organized InAs/GaAs or In- 
GaAs/GaAs quantum dots strain fields result in strong 
piezoelectric fields [9, 10], which are in part responsi- 
ble for the large observed fine-structure splittings, and 
partially mix the deformation potential and piezoelec- 
tric couplings. Although the strength (and even sign) 
of these internal fields is still a matter of intense debate 
[9, 10, 24] and the calculation of the coupling asymme- 
try beyond the scope of this paper, we shall make the 
reasonable assumption that the phonon-exciton matrix 
elements differ by a small phase factor r\ viz. g\ = e lT) g\. 
We shall show next that even small asymmetries 77 can 
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result in a significant entanglement loss. In doing so, we 
exploit the fact that the independent Boson model can 
be solved analytically [11]. Since no transitions between 
different excitonic states are induced by the phonon cou- 
pling (3), the time evolution operator is diagonal [23, 25] 



U(t) 



-Wot 
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( 4 ) 

with a\ = g\/uj\ and D\ the usual bosonic displace- 
ment operator. The time evolution described by eq. (4) 
accounts for the polaron buildup. Let us consider first 
a four- wave mixing experiment [12], where the quantum 
dot is initially brought into a superposition state |0) + |1) 
and the polarization decay |(l|t/(t)/9oe lff °*|0)| through 
phonon dephasing is subsequently monitored. The gray 
lines in the inset of Fig. 2 show results of simulations with 
typical parameters for small quantum dots. The polariza- 
tion amplitude drops on a picosecond timescale to a value 
that is several ten percent smaller than its initial value. 
At later times it remains constant, indicating the forma- 
tion of a stable polaron. It was shown in Refs. [23, 26] 
that this initial drop is due to the emission of a phonon 
wavepacket away from the dot and the resulting imprint 
of the superposition state into the solid-state environ- 
ment. The gray lines in the main panel of Fig. 2 report re- 
sults of simulations where the system is initially brought 
into the superposition state |1) + |2), corresponding to 
the buildup of a photon-exciton entanglement through 
the biexciton decay, see Eq. (1). The polarization am- 
plitude |(2|J7(t)po?7t(£)|l)| slightly drops at early times 
due to the asymmetric phonon coupling through which 
the environment acquires information about the inter- 
mediate exciton state (see panels Xh and Xy in Fig. 1). 
However, because of the small rj value used in our sim- 
ulations the resulting loss of entanglement is extremely 
small (less than a percent). 

Things change considerably if additional phonon cou- 
plings are considered. Indeed, it is well known that 
the polarization amplitude further drops at later times 
[12]. This has been attributed to higher-order phonon 
processes [27], extrinsic charging effects [28], or elastic 
phonon scatterings at the device boundaries [16]. The 
latter are particularly important for single quantum-dot 
devices, as the one considered in this work. Typical 
phonon linewidths due to elastic boundary scatterings 
are of the order of a few tens of peV [16], where in this 
work we shall discard all directional and wavevector de- 
pendencies and use for simplicity a constant value jd- 
We model the elastic phonon scatterings within a master 
equation approach [29] 



-i[H,p]- ld J2[Nx,[Nx,p}}, 
x 



(5) 



with N\ — a x a x - In accordance with elastic scatter- 
ings, this form conserves the number of phonon excita- 
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FIG. 3: (color online) (a) Entanglement decay rate [see 
Eq. (2)] as a function of asymmetry parameter r\ for different 
elastic phonon scattering rates 7<j. The lattice temperature is 
10 K (dark lines) and 20 K (gray lines), (b) Concurrence as a 
function of r/ and for different values of jd and lattice temper- 
ature, and for a radiative decay rate T = 1 ns _1 . The inset 
reports a typical two-photon density matrix of the biexciton 
cascade, where the off-diagonal elements decay according to 
,T . The overall concurrence 1/(1 +7/T) results from this 
density matrix through time averaging over the whole decay 
process. 



tions and gives in thermal equilibrium an exponentially 
damped correlation function (a x (t)a x (Q)) = e^ UJx ~" <d ^ t h\, 
where n\ is the Bose-Einstein distribution. We solve the 
master equation (5) numerically by discretizing the time 
domain into subintervals of length St and expanding p 
and U (St) in the basis of phonon number states. The 
dark lines in Fig. 2 report results of simulations includ- 
ing elastic phonon scatterings, showing a constant decay 
of the polarization amplitude at later time. This is at- 
tributed to the combined effect of polaron formation and 
elastic phonon scatterings. In the formation of a po- 
laron, the exciton becomes partially entangled with the 
phonon degrees of freedom. In turn, the elastic boundary 
scatterings of the phonons contributing to the polaron 
lead to an effective "measurement" of the environment 
on the superposition properties of the intermediate exci- 
ton state. Thereby the system becomes einselected from 
a maximally entangled to a mixed state. Thus, after a 
given time r, when the second photon is emitted in the 
exciton decay, the reduced degree of photon-exciton en- 
tanglement is swapped to a reduced photon-photon en- 
tanglement. 
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The decrease of the polarization amplitude with time 
turns out to be reasonably well described by an expo- 
nential decay with rate 7 [see also Eq. (2)], which is 
reported in Fig. 3(a) for different asymmetry param- 
eters and phonon broadenings jd- One observes that 
even small asymmetries r\ in the phonon coupling result 
in a significant entanglement loss. For instance, using 
7(77 = 0.1) ~ 0.25 ns _1 at 10 K for a moderate phonon 
broadening -fd = 25 (J,eV, we obtain for a typical radiative 
decay rate T ~ 1 ns _1 a decrease of the concurrence by 
about twenty percent [see Fig. 3(b)]. Increased values of 
77, -fd, or temperature even result in higher losses. 

The mechanism described above might explain the 
small degree of photon entanglement observed in recent 
experiments [3] , and is at the same time compatible with 
the long measured spin scattering [30] and dephasing [31] 
times observed in large quantum dot samples. Theoreti- 
cally, the magnitude of the phonon asymmetry might be 
estimated from calculations for realistic quantum dots 
[9, 10], including strain and piezoelectric fields. Experi- 
mentally, a time-resolved measurement of (r) would 
allow a direct observation of concurrence decay. Quite 
generally, the degree of photon entanglement can be en- 
hanced by decreasing the radiative decay times through 
stronger microcavity couplings [19, 32]. In addition there 
might be other coupling channels to the environment, 
such as optical phonons or localized vibrations decaying 
through anharmonic processes, or impurity defects and 
charging centers in the vicinity of the dot [28] that pro- 
duce fluctuating electric fields, which couple differently 
to the exciton states and thereby result in entanglement 
losses. 

In summary, we have theoretically investigated the ef- 
fects of environment couplings in the cascade decay of a 
biexciton confined in a single quantum dot. We have 
shown that a slightly polarization-dependent exciton- 
phonon coupling results, through the formation of a po- 
laron in the vicinity of the dot, in a small loss of photon 
entanglement. This effect becomes drastically enhanced 
when the phonons contributing to the polaron scatter at 
the device boundaries. In consequence, the concurrence 
may easily drop by several ten percent even for moderate 
clastic phonon scatterings and low temperature. Our re- 
sults suggest that in addition to the fine structure split- 
ting, which can be controlled by external magnetic or 
electric fields, also extrinsic effects of the environment 
might provide a limiting factor in quantum-dot based 
entangled-photon sources. 
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